Antiferromagnetic noise correlations in optical lattices 
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We analyze how noise correlations probed by time-of-flight (TOF) experiments reveal antiferro- 
magnetic (AF) correlations of fermionic atoms in two-dimensional (2D) and three-dimensional (3D) 
optical lattices. Combining analytical and quantum Monte Carlo (QMC) calculations using experi- 
mentally realistic parameters, we show that AF correlations can be detected for temperatures above 
and below the critical temperature for AF ordering. It is demonstrated that spin-resolved noise 
correlations yield important information about the spin ordering. Finally, we show how to extract 
the spin correlation length and the related critical exponent of the AF transition from the noise. 
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Atoms in optical lattices hold the potential to unravel 
the fundamental physics of phenomena related to quan- 
tum systems in periodic potentials including spin phases 
and high-T c superconductors One has already ob- 
served the Mott insulator transition for bosons [2|, the 
emergence of superexchange interactions the transi- 
tion between metallic, band-insulator and Mott phases 
for fcrmions 0, @] , and fermionic pairing Q . In addition 
to creating these lattice systems at sufficiently low tem- 
peratures T, a major challenge is how to detect the var- 
ious phases predicted theoretically. These phases can be 
investigated via higher order correlation functions which 
only show up as quantum noise in most experiments. 
Quantum spin noise spectroscopy 0, Q and the measure- 
ment of correlations of the momentum distribution of the 
atoms after release from the lattice (TOF experiments) 
are two ways to probe these correlation functions ||. 
TOF experiments have already been used to detect pair- 
ing correlations in a Fermi gas [l(| , bosonic bunching and 
fermionic anti-bunching of atoms in optical lattices 111 ], 
the Mott-superfluid transition for bosons [T^], and the 
effects of disorder in the Mott phase [l3| . 

Here, we show how AF correlations of fermionic atoms 
in optical lattices give rise to distinct measurable sig- 
nals in TOF experiments even above the critical temper- 
ature for magnetic ordering. Spin-resolved experiments 
are demonstrated to yield additional information which 
can be used to identify the magnetic ordering and bro- 
ken symmetry axis. The main results are illustrated in 
Figs. [l][2] which show noise correlations in the momentum 
distributions after expansion as a function of temperature 
and momentum. We finally discuss how the spin corre- 
lation length and the related critical exponent v can be 
extracted experimentally from the noise. 



We consider a two-component Fermi gas in an optical 
lattice of size N = N x N y N z . In the limit of strong repul- 
sion the gas is in the Mott phase for low T at half-filling 



and can be described by the Hciscnbcrg model 

H=JY. ^™ + S l § rn + (1 + A)«]. (1) 

(l,m) 

Here S; is the spin-1/2 operator for atoms at site r; and 
(l,m) denotes neighboring pairs. The interaction is J = 
At-^ti/U , with U > the on-site repulsion between atoms 
and t a the spin-dependent tunneling matrix element. 
The anisotropy parameter is A = — (tj — ti) 2 /(t 2 - + tV). 
Below we consider both cubic (3D) and square (2D) lat- 
tices with lattice constant a of unity. We do not include 
any effects of a trapping potential. 

A major experimental goal is to detect the onset of AF 
correlations with decreasing T.TOF experiments probe 
correlation functions of the form 111] 
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the atomic den- 
i(r) = 



where A(r) and B(r) are atomic observables measured 
at r after expansion. In Refs. [Tot 
sity correlations C nn (d) were measured with 
^( r ) = S<r^a-( r ) the density operator of the atoms. 
Here n a (r) = ijjy (r)tp a (r) with ip a (r) the field opera- 
tor for atoms in spin state a. We also consider the 
correlations between images of each of the spin compo- 
nents C||(d) corresponding to A(r) = n T (r) and B(r) = 
rij(r). This may be achieved by, e.g., spatially sepa- 
rating the two atomic spin states using a Stern-Gcrlach 
technique [bj]. By applying a 7r/2 pulse before the ex- 
pansion one can also gain access to the spin noise per- 
pendicular to the z component Cj_(d) corresponding 
to i(r) = exp(is y 7r/2)n T (r)exp(-iF7r/2) and B(r) = 
exp(is y TT /2)ni(r) exp(— is y Tr/2). 

After normal ordering and expansion of the field oper- 



2 



ators in the lowest band Wannier states, we obtain 



C nn (k) 

C||(k) 
Gx(k) 
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2^ ~ <W " <*P-k> " <^k>, (4) 
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where K is a reciprocal lattice vector, and = 
A -1 J^i Sje - ' . We assume free expansion of the atoms 
for a duration <, neglect autocorrelation terms oc 5(d) 
in (J3J)- (EJ) : and express C in terms of the momentum 
k = md/t (h = 1 throughout). These correlation func- 
tions have contributions with different scalings. For 3D 
systems, the spins order below the Nccl tcmpcraurc Tjv- 
Assuming the broken symmetry axis along the z axis, 
C nn and C± have a contribution from (s k s k ) ~ 0(1) 
at k = (7r,7r,7r) for T < Tjy- At other momenta or for 
T > Tjy the correlation functions scale as 1/N . Note that 
we will perform calculations where the total spin oc §q of 
the lattice is allowed to fluctuate, which is different from 
the experimental situation where the number of particles 
in each spin state is fixed. This may give rise to different 
momentum independent terms in ([3])- (T5J) for the exper- 
iment, but the k dependent part should be accurately 
captured by our calculations. 

Experiments take 2D column density images of the ex- 
panding cloud corresponding to integrating over z and z' 
in both the numerator and denominator in lp]). The cam- 
eras also introduce a smoothening in the x^-plane, which 
can be modeled by convolution with a Gaussian . In 
total, the experimental procedure corresponds to mea- 
suring the averaged correlation function 



C cxp (k ± ) 
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^T^ {k,± - k±)2 C(k'), (6) 



where k^ = (k x ,k y ) and k = w/l with I = 2irt/ma 
(keeping a for clarity) and w a width depending on the 
pixel resolution of the CCD camera. The averaging in 
the z-direction reduces the contributions to C nn and C± 
at k = (7r,7r,7r) from 0(1) for T < T N to l/N z and the 
Gaussian smoothening further reduces it to 0(1/ k 2 N). 
This reduction happens because the correlations are re- 
stricted to a single point k = (tt, tt, tt) and we are averag- 
ing over N z ■ N x N v /k 2 points. Away from kj^ = (tt,tt) 
or for T > T/v the correlations have a wider distribution 
and they are less affected by the averaging. 

We plot in Fig. Q] 0°* p (fc,fc) as a function of k = 
k x = k y and T for a 3D lattice of size N = 32 3 . The 
spin correlation functions (§k • s_k) were calculated us- 
ing QMC simulations using the Stochastic Series Expan- 
sion method [l5| with directed- loop updates [lj|. This 
method is very efficient for Heisenberg models and gives 
accurate results for a wide range of T for large systems. 
For the isotropic Heisenberg model, the QMC calcula- 
tions yield T N ~ 0.945 J in agreement with Ref.H C°l P 
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FIG. 1: (Color online) Plot of Cnn P (k±) versus temperature 
T and momentum along a diagonal cut k — k x = k y . 



was then calculated from with the k z average in (JSJ) 
included to simulate the experimentally relevant situa- 
Thc main feature of the plot is the dip at k = tt 
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tion 

coming from AF ordering. For T < Tjy , this gives rise to 
a large Bragg dip at k x = k y = tt 0, [TtJ ■ We see that 
the Bragg dip remains also above XV but has a larger 
width due to AF correlations without long range order. 
The spin correlation length £ can be extracted directly 
from the width of the dip which decreases with decreas- 
ing T — T/v reflecting that £ increases as £ ~ 1 / 1 T — Tjy | v 
for T — ► Tjy. We show below how the noise indeed can 
be used to extract the critical exponent v. At high T the 
AF correlations lead to singlet formation of neighboring 
spins. This gives rise to a precursor of the Bragg peak 
at k x — k y = tt and an equal signal of opposite sign at 
k^ = [see ([7])]. These momentum correlations can be 
understood by noting that two fermions in a singlet are 
more (less) likely to have the same (opposite) momentum 
due to the Pauli exclusion principle. The uniform spin 
noise case k = was also considered in Q . 

In the high-T limit J/T = 1/f < 1 [k B = 1), con- 
trolled analytical results for the correlation functions (J3j)— 
([5]) may be obtained by expanding in T~ l . We obtain 
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with Z = 4(6) for 2D(3D) lattices. For simplicity we 
have taken A = and 7k = Z _1 ^ a 2cos(k • a) where 
^2 a sums over the unit vectors spanning the lattice. The 
average of (|7|) and ([8]) over z can be obtained by simply 
omitting the z direction in the sum. 

To obtain analytic expressions for low T, we perform a 
spin-wave calculation, which should be fairly accurate for 
a 3D system at T T^. This yields after some algebra 
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FIG. 2: (Color online) Noise correlation functions C ± xp (kx) 
(red), q xp (k ± ) (black), C£ p (k ± ) (green), and CS p (kx) 
(blue) versus T at kx = (a) and kx = (7r/4, tt/4) (b). 

The dotted (solid) lines are with (without) Gaussian smearing 
(k = 1/40). Dashed lines show the analytical results ©-I©. 
Inset (a): same plot but zoomed in near Tjy. 



where tanh20k = 7k/(l + A) and /k = [expwk/T 1 — 
is the Bosc distribution function for the spin-waves with 

We also have \o ; 



energy w k 



3JV(1 + A) 2 - 7 £. We also have (sgsg) = 
£ k ' sinh- 2 (/3w k /2)/2W 2 where the sum extends over the 
AF reduced Brillouin zone A similar but somewhat 
more lengthy expression can be derived for C±_ (k) . 

We plot in Fig.ElC^kx), C° xp (k ± ), and p (k ± ) as 
a function of T for kx = (7r, it) (a) and kx = (tt/4, 7r/4) 
(b) for the same parameters as above. The solid curves 
include the average over k z only and the dotted curves 
include a Gaussian smoothing in the (k x ,k y ) plane as 
well. The values for T — > with Gaussian smearing are 
consistent with the results of Refs. [ll| given the different 
system size and slightly different measured quantities. 
For simplicity we have excluded the Gaussian averaging 
for kj_ = (tt/4, 7t/4) since it [contrary to the kx = (tt, 7t) 
case] leads to only negligible changes. We trigger the AF 
ordering along the z-direction for T < T/v by including a 
small anisotropy A = 0.01 for C|| xp and C^ xp . 

Figure El^a) shows the Bragg dip for T < Tv coming 

from (s k sf_ k ) ~ l( s f)| 2 f° r ^ = ( 7r ' 7r 7 7r ) Wltn l( s f)l > 0- 
Extrapolating to T = we have \{sf)\ ~ 0.43 ± 0.01, a 
value reduced from 1/2 by quantum fluctuations. On the 
other hand the spin noise parallel to the broken symmetry 
axis C || is reduced with the onset of magnetic ordering for 
T < Tv as can be seen from the inset in Fig. [2ja). (The 
minimal value of C\\ at Tv is dependent on the anisotropy 
A.) Note that even though the density correlations are 



larger and hence more easily measurable, the spin re- 
solved measurements are crucial in verifying e.g. whether 
the correlation dip is indeed due to magnetic ordering 
and not caused, for example, by formation of a period- 
doubled charge-density wave. The difference between Cm 
and C± for A > can furthermore be used to identify 
the broken symmetry axis. In the isotropic case, there is 
no broken symmetry axis and C|| xp = C^ xp = C° s xp for 
A — > 0, as can be seen from the green curve in Fig.[2ja). 

The calculated correlations are rather small even for 
the quantities including the macroscopic contribution 
(*k s -k) ~ l( s f)| 2 at k = (tt, 7r,7r), e.g., for low tem- 
peratures T < J/2 C£P « 10 -2 (10" 3 ) without (with) 
Gaussian smoothening and scales like \/N z (1/JV). This 
is, however, still significantly larger than the correla- 
tions of order C, m ~ 10 -4 which were recently measured 
with slightly bigger system sizes [ljj . The correlations 
are also close to the measured experimental values even 
for temperatures above T/v, i.e., for T = 2 J we have 
Cnn P ~ 6 • 10 -5 (scaling as 1/N), which is comparable to 
the recent experiments fill ]. The measurement of noise 
correlations can thus be used to show the onset of AF 
order even above the critical temperature. Furthermore, 
for current experiments the spin temperature is uncer- 
tain because there are no sensitive probes. If the noise 
correlations are measured, the detailed theoretical curves 
presented here would provide a means of assessing the 
spin temperature in the experiments. 

The density and spin noise at kx = (tt/4, tt/4) de- 
picted in Fig. [2Jb) exhibit a different behavior from that 
at kx = (7T, 7t): it now decreases in numerical value for 
high T and even changes sign above TV- This unusual 
behavior is a geometric effect of the lattice. Note that 
since the noise scales as 1/N away from the Bragg point, 
a measurement requires higher experimental resolution 
than what is presently available or more sensitive detec- 
tion methods such as spin noise spectroscopy [1] . 

Atomic gases are well suited to study fundamen- 
tal problems in 2D physics as the observation of the 
Berezinskii-Kosterlitz-Thouless transition illustrates [l|| . 
Recently, single layer 2D atomic gases have been pro- 
duced which avoids the averaging over z discussed 
above [2(|. We now study the Mott phase at half-filling 
for a 2D square lattice. In 2D there is no ordered phase 
for T > due to fluctuations, but there are still sig- 
nificant AF correlations [2l]]. This is illustrated in Fig. 
EJa), which shows C; e s xp (k_L) and C°* p (k±) as a function 
of T. The AF correlations give rise to a large dip for 
kx = (7r, 7r) both in the density and spin noise which is a 
precursor of the AF ordered state at T = 0. Since there 
is no averaging over the z-direction, the correlations are 
much stronger than for a 3D system above TV. Figure 
[3] illustrates that the noise exhibits the same non-trivial 
features as a function of k and T as the 3D case. 

Critical exponents characterizing continuous phase 
transitions are often difficult to measure. We now demon- 
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FIG. 3: (Color online) (a) Noise correlation functions 
Ci C s X o P (kx) (red) and CS p (kj_) (blue) versus T at k ± = (tt, tt) 
for a 2D lattice. The curves are obtained from QMC simu- 
lations on 64 x 64 lattices. The dotted (solid) lines are with 
(without) Gaussian smearing using k = 1/40. Inset: same as 
solid lines in main panel except at kx = (7r/4, tt/4). (b) Same 
as Fig. [TJbut for a 2D system without the k z summation. 



strate how noise measurements can be used to obtain 
the critical exponent v. The correlation length £ can 
be extracted from the width of the AF dip in the 3D 
resolved function C nn (k) at k= (tt, tt, tt). Wc find that 
this yields a critical exponent v s=s 0.70 as expected 
for a 3D Hciscnbcrg model [22]. Importantly, v can 
also be extracted from the experimentally relevant k z - 
summed correlation function C°* p . Figure QJa-c) show 
C°* p (kjJ at three fixed temperatures above Tjv- To ob- 
tain v, we fit C™ p (kx) to a summed lattice propagator of 
the form £ fc= [2(3 - Z a=xy , z ^s(k a - tt))£ 2 + l]" 1 with 
k z = 2irn z /N z . Figure Bid) shows the extracted correla- 
tion length £ for various system sizes. One clearly sees 
the finite-size effects setting in at decreasing T — T^. The 




extracted power-law yields v m 0.70. To obtain a robust 
value for v one needs to probe the noise for temperatures 
where it is somewhat smaller than what has been mea- 
sured to date. This would however enable one to probe 
the critical properties of the AF transition. 

In summary, we performed analytic and numerical 
calculations modeling TOF experiments for repulsive 
fermionic atoms in optical lattices using experimentally 
realistic parameters. This demonstrated that such ex- 
periments are well suited to detect AF correlations both 
below and above the critical temperature. Spin-resolved 
measurements were shown to yield valuable additional 
information and we finally discussed how to extract the 
critical exponent governing the correlation length close 
to the AF transition from the noise. 

We thank H. Moritz for the suggestion to look at TOF 
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wegian Research Council. B. M. A. acknowledges support 
from the Villum Kann Rasmussen foundation and E. D. 
from NSF DMR-0705472, CUA, DARPA, MURI. 



FIG. 4: (Color online) Panels (a-c) show (-l)C„f 1 p (k x ) at 
T/J = 1.0 (a), T/J = 1.2 (b), and T/J = 1.6 (c). (d) 
Log- log plot of extracted correlation length £ versus \T — Tn\ 
displaying the power-law behavior £ ~ \T — Tjv|~ 0,70 . 
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